Real-time renormalization group and cutoff scales in nonequilibrium applied to an 
arbitrary quantum dot in the Coulomb blockade regime 
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We apply the real-time renormalization group (RG) in nonequilibrium to an arbitrary quantum 
dot in the Coulomb blockade regime. Within one-loop RG-equations, we include self-consistently 
the kernel governing the dynamics of the reduced density matrix of the dot. As a result, we find 
that relaxation and dephasing rates generically cut off the RG flow. In addition, we include all 
other cutoff scales defined by temperature, energy excitations, frequency, and voltage. We apply the 
formalism to transport through single molecular magnets, realized by the fully anisotropic Kondo 
model (with three different exchange couplings Jx, Jy, and Jz ) in a magnetic field . We calculate 
the differential conductance as function of bias voltage V and discuss a quantum phase transition 
which can be tuned by changing the sign of JxJyJz via the anisotropy parameters. Finally, we 
calculate the noise S(fl) at finite frequency Q, for the isotropic Kondo model and find that the 
dephasing rate determines the height of the shoulders in dS{yi)/dQ, near Q. = V. 

PACS numbers: 73.63.Nm, 05.10.Cc, 72.10.Bg 
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I. INTRODUCTION 

A fundamental issue of recent interest is the de- 
velopment of renormahzation group (RG) methods in 
nonequihbriumiii2ii2iiii^i^i^i§ Besides the discovery of new 
power law exponents for the conductance induced by 
nonequilibrium occupation probabilities an interesting 
question was raised whether voltage-induced decay rates 
provide additional cutoffs of the RG flow.^'^° In this 
context, the nonequilibrium Kondo model has been dis- 
cussed, which can be realized by a single-level quantum 
dot (QD) in the Coulomb-blockade (CB) regime cou- 
pled via spin exchange processes Jaa' to two reservoirs 
a — L,R. In the isotropic case and above all cutoff 
scales, the exchange couplings J — Jaa' are all the same 
and are enhanced by reducing the band width A of the 
reservoirs according to the poor man's scaling equation 
^ = 2J^," with / = ln(Ao/A) (Aq denotes the inital 
band width). The enhanced screening of the dot spin 
leads to the Kondo effect with unitary conductance below 
the Kondo temperature Tk = Aq exp[— 1/(2J)] (Ref. 
(for experiments in quantum dots see, e.g., Refs. Il3r ). 
However, for voltages V ^ Tk, it was argued that the 
system cannot reach the strong coupling fixed point since 
the nondiagonal coupling constants, Jlr = Jrl, are cut 
off by the voltage"^^ and the diagonal ones, Jll and Jrr, 
by the voltage-induced decay rate F = TTjfj^\A=vy'^ 
This has raised the fundamental question how decay pro- 
cesses can be implemented in nonequilibrium RG. Apply- 
ing flow equation methods to the isotropic Kondo model 
without magnetic field, it was shown within a two-loop 
formalism in Rcf. 3 that the inclusion of a third-order 
term ~ J'^ in the RG equation leads to a cutoff of the 
RG flow at the scale F. 

In this paper, we will analyze this problem from a 
more general point of view and will show within a mi- 
croscopic one-loop RG formalism that relaxation and de- 



phasing rates will always cut off the RG flow for an ar- 
bitrary quantum dot in the Coulomb blockade regime. 
This confirms the conjecture of Refs. andlTo! and gen- 
eralizes the analysis of Ref. 7 to an arbitrary QD (in- 
cluding orbital and spin fluctuations, many levels, inter- 
ference effects, etc.). We propose to use the real-time 
RG (RTRG) formalism of Ref. [U with a cutoff defined 
in frequency space since this approach directly discusses 
the time evolution of the reduced density matrix of the 
dot via a kinetic equation. Within this formalism, the 
decay rates occur naturally as the negative imaginary 
parts of the eigenvalues of the kernel determining the 
dissipative part of the kinetic equation. This has already 
been demonstrated previously by applying RTRG to the 
calculation of steady-state transport through quantum 
dots in the charge fluctuation regime^ and to the study 
of the real-time evolution of the occupation probabilities 
within the spin boson model."^ Another advantage of the 
RTRG approach is the fact that the kernel can easily be 
inserted self-consistently into the one-loop RG equations 
of the coupling parameters (analogous to self-energy in- 
sertions within Green's function techniques), providing 
the unique possibility to obtain the physical decay rates 
within a nonequilibrium one-loop RG formalism. In addi- 
tion, we also provide a microscopic formalism from which 
all other standard cutoff scales, such as temperature, en- 
ergy excitations (e.g., magnetic fields), frequencies, and 
voltages, can be deduced analytically. 

The original RTRGi was formulated with a cutoff de- 
fined in time space for the reservoir correlation function. 
This makes it technically difficult to apply the formal- 
ism to problems where the interaction between dot and 
reservoirs is nonlinear as it is the case for quantum dots 
in the cotunneling regime, where orbital and spin fluctua- 
tions dominate transport. Therefore, we use in this work 
a cutoff defined in frequency space but adapt the same 
formalism to set up the RG equations as in Ref.fT. This 
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leads to a combined time-frequency formalism since the 
time-ordering of the renormalized vertices is needed due 
to their operator nature (the degrees of freedom of the 
dot arc not integrated out within RTRG, and therefore 
all coupling vertices are operators acting on the dot de- 
gree of freedom) . The only disadvantage of the analytic 
formalism presented in this work is still the fact that the 
irrelevant prefactors of the various decay rates cutting off 
the RG flow cannot be determined unambigiously; this 
has to be left for future developments. 



First, we apply the formalism to quantum transport 
through single molecular magnets (SMM). Recently, it 
has been shown that the study of Kondo physics can 
be used for transport spectroscopy of SMMyi^ii^ i.e., the 
various anisotropy parameters determining the spin ex- 
citation spectrum can be identified. In the regime where 
the Kondo temperature is smaller than the distance to 
the next spin excitation, it has been shown that a pseudo- 
spin- 1/2 model can be derived which can be mapped onto 
the fully anisotropic Kondo model with three different 
exchange couplings J^, Jy, and J^. Interestingly, this 
model reveals a quantum phase transition by changing 
the sign of JxJyJz, separating the flow to the weak and 
strong coupling regimes. Since the exchange couplings 
depend on the transverse anisotropy parameters, which 
in turn depend on the coupling of the SMM to the leads, 
this phase transition can be tuned in an experimental 
setup. Using RTRG, we calculate the differential conduc- 
tance G{V) as function of bias voltage at finite magnetic 
field and show that the Kondo-enhanced conductance at 
V = h, where h is the level spacing between the ground 
state and the first excited state, disappears by tuning the 
system through the phase transition. 



Second, we calculate the quantum noise S{fl) as func- 
tion of frequency for the isotropic Kondo model at finite 
bias and zero magnetic field. We find that the dephas- 
ing rate can be identified by studying the derivative of 
the noise near Q ~ V. Specifically, it turns out that 
the noise has a dip at V = (see also Ref. T^, where 
the noise has been calculated for the Toulouse point), 
whereas the derivative dS{^l)/dft shows a characteristic 
shoulder with a height depending on the dephasing rate. 



The paper is organized as follows: In Sec. |TT1 we set up 
the general model and show the relation to the nonequi- 
librium Kondo model. Section IIIII summarizes the dia- 
grammatic language in Liouville space. Section [IV] is the 
central technical part where we set up the RG equations 
and explain how decay rates cut off the RG ffow. Finally, 
we apply the formalism in Sec. |V] to transport through 
single molecular magnets and in Sec. IVII to the calcula- 
tion of quantum noise. Two appendices provide further 
details of the RG formalism. 



II. MODEL 

We consider an arbitrary quantum dot coupled to 
reservoirs via tunneling processes, 

H = H,,s +Hu + Ht. (1) 

where iJics, ^^d, and Ht denote the Hamiltonians of the 
reservoirs, the dot, and the tunneling, respectively. 

H^cs — 'y ' -^rcs ~ 'y ] ^ka<j(^\a^O,kacr (2) 

describes the noninteracting Hamiltonian of the reser- 
voirs with a\^^ (akaa) the creation (annihilation) oper- 
ators, a is the reservoir index, a denotes the spin, and k 
is an index for the single-particle states in the reservoirs. 
Each reservoir is assumed to be infinitely large and de- 
scribed by a grand canonical distribution with electro- 
chemical potential iia and temperature T. The isolated 
dot Hamiltonian is written in diagonalized form as 

Hu = J2e,\s){s\, (3) 

s 

where s is an index for the many-body eigenstates of the 
dot with energy eigenvalues Es. Finally, the interaction 
between dot and reservoirs is described by the standard 
tunneling Hamiltonian 

Ht^Y. ^kiaL.ci. + H.c, (4) 

akl{T 

where c/o- annihilates a particle with spin a in the single- 
particle level I on the dot and t^f denotes the tunneling 
matrix element. 

Since the reservoirs are infinitely large, we describe 
their spectrum by the continuum density of states 
p^(uj) = X)fe '^(^ ^ + A^a); with /J, = aa an index con- 
taining the reservoir and the spin index (this will be used 
implicitly in the following). For the general discussion, 
we include the case of spin- and frequency-dependent 
density of states in the reservoirs. We introduce the con- 
tinuum fields 

a^,+ {uJ)^—^^=y^5{u-ek^, + ^io,)al (5) 
VPm(^) k 

and a^_(a;) — ap_|_(u;)T which fulfill the anticommutation 
relation {a^,,(a;), a^^/,,' (w')} = Sr,-,,i5fj,fj,i5{uj - uj'). With 
this notation, the reservoir Hamiltonian and the tunnel- 
ing part can be written as 

Hrcs = ^ duj {uj + Ha)a,_t+{uj)a^^{u}), (6) 

Ht ^ Y duja^+{u})g^{uj) + H.c, (7) 

with 

5m H = \/Pt^i^) '^tfi^yia, (8) 

I 
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where t""'{Lo) = f^f is the tunnehng matrix element in 
the continuum notation evaluated for reservoir state k 
such that ekaa^fJ'a = ^- The contraction of two reservoir 
field operators with respect to the equilibrium reservoir 
distribution is given by 

(o^^(w)a^'^/(w'))rcs = Sn-,,'6^^>6iuj - uj')d^f^, (9) 

where /+ = and / - = 1 - ^ /_^, with = 
l/[exp(/3ci;) + 1] denoting the Fermi function (note that 
the different electrochemical potentials of the reservoirs 
occur in our notation via the interaction picture from 
the time-dependence of the field operators). 9^ = 9{Ao — 
Lu) contains the initial band width Aq of the reservoirs 
(which are assumed to be all the same relative to the 
corresponding electrochemical potentials). 

We now consider a quantum dot in the Coulomb block- 
ade regime, i.e., the total charge is fixed, and only cotun- 
neling processes via virtual intermediate states can lead 
to orbital and spin fluctuations. The effective Hamil- 
tonian in this regime is standardly derived using the 
Schrieffer- Wolff transformation,^^ leading to 



H, 



off 



H„ 



off, 



(10) 



with 



- 9^^,'i^^^')<^f^'-i^')a,,+ {(^)] , (11) 



where 
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UJ + Ha + Es - Hu Uj' + fia' + Es' - i?D 



(12) 



corresponds to virtual processes, where the electron first 
hops from the reservoir to the dot and then back, and 
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1 



1 



U} + Ha ~ Es 



Hr 



^' + Ma' - Es+ Hy:, 



(13) 
3m('^)|s') 



describes the reverse process. Processes where two elec- 
trons hop on or off the dot are not written here but 
can easily be incorporated (they are only important for 
molecular systems with negative Coulomb interaction; 
see Ref.Hl). 

In normal-ordered form (with respect to the equi- 
librium reservoir distribution), denoted by the symbol 
we get from Eqs HU]), ([11]), and ® 



Eire 



eflf 



off 



(14) 



with a renormalized dot Hamiltonian 



rrcn tt 



duje^T^gl(u,uj)fZ, (15) 



which can contain logarithmic energy renormalizations of 
the dot states due to orbital interferences or due to spin- 
dependent tunneling matrix elements (see, e.g., Refs. [l^ 
and '2(f) (for the Kondo model under consideration in this 
work, such renormalizations do not occur). Finally, the 
normal-ordered interaction term reads 



V =.V^s- = E / dujduj' 
MM' J-^o 

g^^,{uj,uj') : a^+(w)a^'_(tj') :, (16) 



with 



g^^,{uj,uj') ^^gl^,{uj,uj'). 



(17) 



Equations p4|) and (fT6|) are the final general form 
of the model under consideration, which is the starting 
point for the renormalization group formalism. It is still 
completely general, except for the fact that the dot is 
assumed to be at fixed charge. In order to simplify the 
notation, we omit in the following the index "eff", and 
use the short-hand notation 



V = giv : ai+ay 



(18) 



where we sum/integrate implicitly over the indices 1 
wi/ii and 1' = wi'Mi'. We note the property 



511' = ffi'i> 



(19) 



which guarantees the Hermiticity of i?D- 

The fully anisotropic Kondo model under considera- 
tion in Sec. |V] is realized for the special case where the 
dot Hamiltonian consists of two states with (pseudo-) 
spin up or down (i.e., s = ± =t, i denotes the dot spin). 
The antiferromagnetic exchange processes between the 
dot spin S_ and the reservoir spins are described by the 
coupling 



511' 



E 

i—x,y,z 



J I Ql I 



(20) 



where cr*, i — x,y,z, are the Pauli matrices. Insert- 
ing this into Eq. pB]) gives the standard form of the 
anisotropic Kondo model. 



V 



s=± 



(21) 
(22) 



(cj)ap'_(cj') 



where h denotes the effective magnetic field in the z direc- 
tion. In contrast to the usual case — J^, we discuss 
here the fully anisotropic Kondo model with three dif- 
ferent exchange couplings, a model of recent interest if 
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the dot is replaced by a single molecular magnet .J^ In 
this case, the isolated molecule is described by the spin 
Hamiltonian 

n 

(23) 

where D and B2n denote the longitudinal and trans- 
verse anisotropy constants, and the original spin Sm is 
greater than 1/2. hz denotes the physical magnetic field 
in the z direction. If one projects an isotropic exchange 
{J /2)S_^.jg_^^,ay^+{uj)a^i -(uj') between the molecule and 
reservoirs onto the two lowest eigenstates |±) of iJ,„oi 
(which is justified when the Kondo temperature is lower 
than the first magnetic excitation), one obtains a pseudo- 
spin-1/2 model described by the fully anisotropic Kondo 
model [Eqs. ^ and (HH)] with 

h = (+|if,„oi|+) - (-l-ffmoih), (24) 
J^'y = J{+\S+±S^\-), (25) 
J- = 2J{+\S^\+) > 0, (26) 

where J is the isotropic exchange constant between the 
original molecular spin and the reservoirs (see Ref.fTi for 
further details). 

III. PERTURBATION SERIES 

We aim at calculating the stationary dot distribu- 
tion p^*", the stationary current in lead 7, and the 
frequency-dependent noise power 

= sit+s^;-^-2nS{n)i:j:i (2?) 

with Sr = r - 4 and 

S'n'' -I f dte-^^mr{t),r'{0)}). (28) 

Due to current conservation, we have 

E^o"'-E^o"'=0, (29) 

7 7' 

and therefore, for two reservoirs, it is sufficient to calcu- 
late the diagonal noise S^''. 

The current operator H for lead 7 is given by — 
— {d/dt)N^ = ~i[V, N^], where N-y is the particle number 
in reservoir 7 (we use units e = h = I). Using Eq. psp . 
this gives 

P = i{S^a^ - Sja[)9ii' ■ ai+ai'- : . (30) 

Following Ref. IT!, we start from an initial distribution 
p{to) = p{to)prcs which factorizes into an arbitrary dot 



part p(to) and an equilibrium grand canonical distribu- 
tion Pros = IIq exp[— /3(iJr"g — fj.aNa)]/Za for the reser- 
voirs. The reduced density matrix of the dot at time t 
can then be written as 

p{t) = Tr,ese-^^(*-*«)p(io)Pros, (31) 

where L — [H, ■] is the Liouville operator, which is a 
superoperator acting on ordinary operators b via Lb — 
[H, b] . According to Eq. (fT4|) , we decompose L ~ Lies + 
Ld + Lv- Using Eq. (|18p . the Liouville operator of the 
interaction part can be written as 

Ly = [U,.]=p'GP^; : Jf+Jf:_ :, (32) 

where we sum implicitly over the Keldysh indices p, p' = 
±. Here, 

^11' ~ Spp'G^^t (33) 

and Jf^ are superoperators acting on usual dot (lead) 
operators b via 

G++b = gn'b, G--b=~bgiy, (34) 

Ji+ 6 = ai^6, J^^b = bai^. (35) 

Taking matrix elements with respect to the dot states, 
the superoperators Ld and G^f/ are given by 

(-^d)sis;,S2s;j = {HYi)sis.2^s\s'^-' 5a^S2{HYi)s'^s{, (36) 
{Gtv)sis'i.S2s'^ — [giv)siS2^s{s'^, (37) 
{Giy)sis{,S2s'.^ = -^sisal^llOs^s'j- (38) 

If the states \s) are the eigenstates of H^) with eigenvalues 
Es, we get 

(-^d)sis'i,S24 = {Es^ - -Bs'J'^sisa'^s'i^- (39) 

From these matrix representations, we get 

s p s 

which is an important property guaranteeingthe conser- 
vation of probability ^gP{t)ss = 1 (see Ref. Hf). 

Following Ref. |l|, we expand Eq. ([5T|) in Ly and inte- 
grate out the leads in order to get an effective descrip- 
tion for the dynamics of the dot. We define the interac- 
tion picture of Ly with respect to Lres + Ld and obtain 
Lv{t)^p'G\i, ^ : Jl+4_ :, with 

G\l,^ = e*("i-"i+^°i-^"'i)*e*^°*Gf^,'e-'^°*. (41) 

Each term in the perturbation expansion is then aver- 
aged over the equilibrium reservoir distribution by using 
Wick's theorem. Using Eq. ([9]), this leads to pair con- 
tractions between the superoperators Jf^ given by 

(42) 
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current in Laplace space 




FIG. 1: Example of a sequence of two irreducible blocks. The 
double line represents time propagation in Liouville space of 
the dot (time increases to the right). The black dots represent 
the vertices G^^, . The lines connecting the vertices are the 
reservoir contractions. Whereas the auxiliary vertical lines a 
and c hit reservoir contractions, line b does not and separates 
the two irreducible blocks. 



In this way, we obtain a sequence of time-ordered dot su- 
peroperators {—i)G^y ^ in interaction picture, connected 
in an arbitrary way by lead contractions (for details 
and diagrammatic representations on the time axis, see 
Ref.E). 

The series of all diagrams can be grouped in irreducible 
and reducible parts, where irreducible means that any 
vertical cut to the time axis hits at least one reservoir 
contraction (see Fig. [1] for an example) . We define the 
kernel Eq = /o°° (it e*^*I](i) in Laplace space, where T,{t) 
is the sum of all irreducible diagrams between time and 
t. The whole series of all diagrams can then be formally 
resummed and we obtain the following result for the dot 
distribution in Laplace space: 



Pn = 'nnp{to), Llf 



Q — Lr, — iTic 



(43) 



with pii = dt e^^^*'p{t) . The stationary distribution 
follows from p^^ — — zlimfj^o ^PQ, leading to 



(LD + iS)p^* -0, 



(44) 



where E = Sq^q- Therefore, the central quantity to be 
calculated within renormalization group is the irreducible 
kernel S; the stationary distribution then follows from 
finding the eigenvector with eigenvalue zero of Ld + 
We note that the irreducible kernel starts and ends with 
two boundary vertices, denoted by B and A, respectively. 
Compared to Eq. (|4ip , their interaction picture is slightly 
differently defined and contains the frequency Q, 



APP 

^ll'f2.t 



n,t 



'^'^ A\\,^e'^^-\ (45) 
iOtg«('^i-'^l+Pai-A'Q,'j)* e*"^°*i?^f'j2 (46) 



Before starting the RG, we have Aff = i^if/n = G^*^, , 
but during RG, the boundary vertices renormalize differ- 
ently and can become 17 dependent. 

A similiar approach can be set up for the calculation 
of current and noise. Choosing to = 0, we write for the 



where T is the time-ordering symbol and 



(47) 



Ll = -^{I\■}^p'G'l^^,^ : Jf+Jf,_: (48) 



is the current superoperator. The current vertex in Li- 
ouville space is given by 



with 



(49) 



(50) 



Expanding the exponential as described above in Ly, 
integrating out the reservoirs, and resumming the whole 
series using the irreducible blocks, one arrives at 



(51) 

and the stationary current follows from 

4l = TrDS>^*, (52) 

with S'' = Sjj^Q. Here, Trn denotes the trace over the 
dot states and is the irreducible kernel containing 

exactly one current vertex Gjy^ with interaction picture 
defined by 

7,PP' _ piatpi("i-'^'i-l-A'ai-/i„'^)* piLDt(^7,Pp'p-iiDt 



r<y'PP 



(53) 

For the noise, we choose to — —oo and start from the 
expression 



C77 



dtf 



TrTHLimHLl(t)] 
g-,:/°^<it'Lv.(t')p(_oo)^^^^. (54) 



Again, expanding in Ly, integrating out the reservoirs, 
and resumming via irreducible blocks gives 

S^-^' =TM^l^' +i:lUn^f)f\ (55) 

with T,J^ the irreducible kernel containing exactly two 



current vertices (see also Ref. l21f ) . Since the current ver- 
tices can also lie at the two boundaries of the kernel, 
one has to define several boundary current vertices with 
slightly different interaction picture compared to Eq. 
(see Appendix B for more details). 
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IV. RENORMALIZATION GROUP 
FORMALISM AND CUTOFF SCALES 

We now take a reduced band width A in the definition 
of the contraction [Eq. ((42|)] by replacing 0^ 9{A— 
Following Ref. [H, we determine the A-dependence oi^ Lb 
and Gil' in such a way that the total sum of all diagrams 
remains invariant. This leads to the RG diagrams of 
Fig. [2] which are evaluated in Appendix A with the result 




FIG. 2: The RG diagrams determining the renormalization of 
(a) G and (b) Lu- The cross indicates differentiation of the 
lead contractions with respect to A. Time-ordering is defined 
with respect to the middle time. 



rPlPl 



dA 



dtS^ 



-P2/5(G?iri/2)..(G?^!':./2).4' (56) 



(57) 



with 5lj = 5{A — {G)ij — {i\G\i) denotes the matrix 
element with respect to the eigenvectors of Ld , 



L 



A, 



(58) 



h and F > describe dot excitations and decay rates, 
respectively. Since F leads to exponential damping be- 
tween the vertices, the time integrals can be cut off by 
F^^, thereby neglecting only small perturbative correc- 
tions for energy scales below F. Note, however, that Ld 
has a unique eigenvector |0) with zero eigenvalue since 
the system is approaching a stationary state.-^ There- 
fore, the contribution from this eigenvector does not lead 
to exponential damping but it will be shown below that 
it does not contribute to the RG flow in leading order. 
For later purpose, we note that the "ket" form (ss' |0) de- 
pends on the specific problem under consideration, but 
the "bra" form (0|ss') is unique and is given by 



1 



(59) 



where Z is the number of many-particle states considered 
on the dot. This property follows directly from Eq. 

Using Eqs. (|4T|) and ((58)) . we obtain the following ex- 
pression for the time integral in Eq. l[56|) [± corresponds 
to the two terms on the right hand side (rhs)]: 



Jo 

(60) 



±(w2 + x±) - i{Tj - Tik) 
with \u}2\ = A, \k = (A^ + Afe)/2 = h^k - i^ik, x± ^ 
and wii' = (wi -I- uj[)/2. This provides a cutoff at 



A± — max(|a;±|, Fj, |Fj — Tik\), containing frequencies, 
voltages, dot excitation energies, and decay rates. Above 
the cutoffs, we obtain ±sign(w2)/A for Eq. ([SO)) . There- 
fore, we can replace — sign(w2)p/^(w2) by sign (w2)p[ 1/2 — 
/P(w2)] = 1/2 - /(A) « (1/2)6*7 in Eq. 1^, providing 
the cutoff set by temperature. This gives 

+ {Gi2^^ )ij {G211 ^ )jk 

-0A_(G^fr'O.,(GfD,fe}, (61) 

where I = ln(Ao/A) denotes the flow paramter. First, 
we get from this equation the central result that decay 
rates always lead to a cutoff of the RG flow. If all k 
correspond to the eigenvector with zero eigenvalue, we 
get A_(. = A_ and the two terms on the rhs of Eq. (pT]) 
cancel. If at least one eigenvector has nonzero eigenvalue, 
we obtain a cutoff either from F^- or \Tj — Fifc|. Neglect- 
ing the irrelevant difference between the various decay 
rates, we replace them in the following by an overall scale 
F. Second, above all cutoff scales, the RG equation pre- 
serves the initial form of the vertex given by Eqs. ([55]) 
and Below, we show that a similiar analysis leads 

to Ld — [Hu, ■] in leading order, with i/o = J2s 
denoting the renormalized dot Hamiltonian [see Eq. (j66p ]. 
Inserting these forms in Eq. ((6T|) . we can project the RG 
equation for the vertex on one part of the Keldysh con- 
tour and we obtain the final result 



{dgiv/dl)ss' = ^2^'^2 

X {^max(T,|2;+|,r) {912) ss{g2v) ss' 

- ^max(T,|K_|,r) {92v)ss{gi2)ss'] , (62) 

with X± = ^Ia2 - Maia'i ^ ^11' ± i^S - E^s') and Ess' = 

{Es + Es>)/2. If the frequency dependence is irrelevant, 
we get gfj,fj,'{LUi,Lu[) = g^^'{ujii'). This equation is a gen- 
eralization of the RG equation of Ref. 5> to an arbitrary 
QD in the CB regime with (possibly) frequency depen- 
dent density of states in the leads, including the micro- 
scopically derived cutoff scales from decay processes. 

We now turn to the leading order analysis for the RG 
equation ([57)1 of the dot Liouvillian. We first insert the 
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leading order form (j33p for the vertex and get 

Using {d/ dK)6^^9uj2 — OujiSuj2 +Soji(^uj2 ^^'^ interchanging 
1 ^ 2 in the second term, we get 




+ f^rf^AGZ/2hiG72.-t/2hk}- (64) 

Performing the same steps as for the derivation of 
Eq. (pT|) . we obtain in leading order 

x{0A+p7^;(G?f).,(Gf'/U + 

+ 0A_p^7(G^;p'),,(G?f),,}, (65) 

with A± = max(|x± I, Fj, jFj — Tik\) and x± — ^1^2 ~ 
— wi ± {hj — hik). Analogous to the conclusion drawn 
from Eq. ((6T|) . we see that decay rates will always lead 
to a cutoff of the RG flow for L-q. Here, the case that 
1 = corresponds to the eigenvector with eigenvalue zero 
can be excluded, since X]p(0|Gi2 = due to Eq. ([55]) 
and the property (j40p which is conserved under the RG 
flow. Second, due to the leading order form ([M)) of the 
vertex, G^^ = gi2- acts only on the upper part of the 
Keldysh contour and Gj~2~ = — ■ 512 only on the lower 
one. Therefore, they commute, and above all cutoff scales 
(i.e., for A ^ A±), we obtain no contribution from p' = 
—p to the renormalization of Ld . The contribution from 
p' — p gives the leading order form L^' = [Hy) , •] with a 
Hcrmitian renormalized dot Hamiltonian Hj^ . In analogy 
to Eq. (p^ . the RG equation for Hd reads 

{dHu/dA)ss' = 

X {^max(T,|x+|,r)/i (5l2)ss (ff2l)gs' 

+ ^max(T,|£c_|,r)/Ji {92l) ss{gi2) ss'} , (66) 

with x± = /iQ2 — /ictj — wi ± {Eg — Ess')- If the frequency 
dependence of 512 is irrelevant, we obtain in leading order 

{dHu/dA)ss' = 26'inax(T,|y|,r) ig^l^^')ss{9^L'^l)ss', (67) 

with y ^ ~ fJ-a' - Eg + Ess' ■ 

The RG equations ((62)) and ([66)1 are the central results 
of this section. They provide the leading-order renor- 
malization of the vertex and the dot Hamiltonian for a 
generic quantum dot in the Coulomb blockade regime. 



Besides the full frequency dependence and the influence 
of temperature and voltage, they include the influence of 
the renormalized dot energies and the decay rates on the 
RG of the vertex. The renormalized dot energies follow 
from Eq. ([66]) but we still have to set up the RG equation 
for the decay rates Vi . They follow from Eq. , where 
we insert on the rhs the leading order form ([55]) and ([M]) 
for the vertex, and the leading order form ' = [Hd , ■] 
for the dot Liouvillian, i.e., we neglect essentially the in- 
fluence of the decay rates on themselves. As shown above, 
the leading order form L^' = [^^d, ■] arises from the prin- 
cipal value part of the time integral and taking p — p' . 
There are two additional contributions to Ld. The first 
one arises from p — p' but taking the (S-function part of 
the time integral. This leads to a contribution of the form 
= {i/p' , •} with an anti- Hcrmitian dot Hamiltonian 
describing energy broadening. The second one arises 
from p = — p', i.e., from diagrams connecting the upper 
with the lower part of the Keldysh contour. This part 
is denoted by ^"^^ describes the physics of relaxation 
and dephasing. For p = — p', only the (^-function part of 
the time integral contributes and a straightforward cal- 
culation gives the results 

- i ( -TT- ) = -271- 9t,t,' (w, UJ')s^s29^i^.' i^J, ^')\s'2 

y-^{e^e^')f-f+ 5{lo-u' + y) (68) 

together with 

{H^^)ss'^{-l/2)Y,{U^)-s-s,s's, (69) 

s 

where y — ^J■a — fia' + Ei ~ E2 and Ei — Es-s'.- The 
frequency integrals over uj and lu' can be calculated an- 
alytically due to the two (5-functions. If we take the 
frequency dependence of the vertex in leading order 
g^^'[uj,uj') = 5^^'(^^^±^), we obtain explicitly 

V / Sis'^,S2s'2 pp' = ± 

X5PM'(P'(A-I||))s,s. wb'(A-|||)):i4. (70) 

Thus, for T = 0, we see that decay rates are only gener- 
ated for \y/2\ < A < \y\, i.e., essentially below all cutoff 
scales. The reason for this interval is a simple golden rule 
argument illustrated in Fig. [3] for the case s = Si = s[, 
s' = S2 — s'2 and Es = Es' ■ Finally, we note that includ- 
ing the influence of the decay rates in the rhs of Eq. (|68p . 
one obtains essentially a broadening of the (5-function by 
F. 

A similiar analysis can be used to determine E, , 
and YP'^ . However, since several boundary vertices with 
different cutoff scales have to be distinguished for the 
general case, we summarize here only those matrix el- 
ements necessary for the Kondo model [Eq. (HH)] (for 
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A 



ML -f 

V 

MR - 



(a) A > "1/ (b) y > A > V/2 (c) V/2 > A 



tdA 



FIG. 3: Generation of a transition rate from state s to s' 
by the RG (for simplicity, we consider two leads L, R with 
V = — Mr S'lid the case Es = E^i here). Energy-conserving 
transitions from one of the intervals dA which are integrated 
out to the other lead or vice versa can only occur in situa- 
tion (b). 



more details, see Appendix B). The RG equation for 
^ss.s's' = Wss' is identical to the rhs of Eq. ([70]) for 
s ^ s' (note that dA < 0, so that the renormalization 

dWs,' > 0), 



pp'=± 

x0\y\/2O{py)fy^^^^. (71) 



If the nondiagonal matrix elements of the stationary 
distribution are zero, the diagonal probabilities pf — pf^ 
follow from the rate equation 



^ (W^s-pf, - Ws-sPf) = 0, 

s' , s'^s 



(72) 



and the stationary current can be written as 

(73) 



The RG for the current rate J^s ^1' = ^ls,s's' is 
given by the rhs of Eq. ([70|) but multiplied with 2c^^, = 
— {S-ya — S-ya') and summing over s. 



E(^),,, = -2-El5..'b'(A"l|l)WP 

5 \ / ss spp' 

^'^cl^'0\yW20{py)fUZM- (74) 



To calculate the noise at finite frequency f2, we need 
the frequency dependent vertices givn, arising from the 
renormalization of the current vertex [Eq. ([55]) ]. They 
follow from Eq. (p^ by replacing the product of the two 
g vertices on the rhs by the average (l/2)((7n5 + 53n) and 



shifting the cutoff by =Ff2/2 for the two terms of this 
average, respectively. 



{dgivn/dl)s 



1 



X {^'max(T,|x+-f2/2|,r) (ffl2o)ss (521' )ss' 

^ ^max(T,|i;_-fl/2|,r) (521'o)s5(.9l2)ss' 
+ ^'max(T,|2:++0/2|,r) {gi2) ss{g21' il) ss' 
- ^max(T,|2;„+0/2|,r) {g2v) ss{gi2n)ss'] ■ (75) 

For the Kondo problem without magnetic field, it can 
then be shown that the diagonal noise follows from 
^li = \ Yl,ss'i^n')si^' ^ where the RG of the noise rate 
Es.'(l^rr)-' - i:ss'i^'S)ss,s's' + -O) is given 

by the rhs of Eq. ([70)1 but multiplied with 2(c^^,)^, sum- 



ming over s, s', replacing g 



^ \ dA 



gn, y ^ y + ^, and adding 
(76) 



ss' pp' 



x9^y+n\/20{p{y + n))flf^ 



-p 

A-|y+0 



The noise contribution from Sj^ can be shown to be ir- 
relevant without magnetic field, since = i/D, in this 
case (see Appendix B). 



V. TRANSPORT THROUGH SINGLE 
MOLECULAR MAGNETS 

We now apply the formalism to transport through sin- 
gle molecular magnets described by the pseudo-spin- i 
dot Hamiltonian (pij) and interaction (P^ . where 5^^' — 
(1/2)J^„,S'V^^,. The effective magnetic field h [differ- 
ing from the physical magnetic field hz used in Eq. (I23p ] 
and the exchange interactions are given by Eqs. 
Using the form of the interaction in Eq. ([66]) . we obtain a 
constant, i.e., no contribution to Ld- Thus, the energies, 
given by Eg — sh/2, s = ±, stay invariant. From the 
vertex RG equation (j62p . we obtain in leading order 

P'^ci^) = 1(0+ + ei){Ji^4^, + JLjL')^ (77) 

where i,j,k are all different, and we have defined 



^max(T,\x±h\X)' 



61 



± — fmax(T,|2;±h/2|,r) 



(78) 



with X = [la — Mqq' ^ ^- From Eq. (j7ip . we obtain for 
the rates (s ^ s') 



dA 



Wss' = -(^/4) E '5L'^i,"',^(roaa' 



x/I/;V/««'b'(A-|^/r<:,|)]^ (79) 



9 



0.06 



0.01 







= 0.006, - 




B2 


= 


0.05 








Bi 


= 0.005, - 




B2 


= 0.005 


0.04 


Bi 


— 005 ■■■ 




B2 


= 0.002 


0.03 


- Bi 


= 0.005, - 




B2 


= 


0.02 













0.5 



1 

V/h 



1.5 



0.8 



0.4 



■oN 



-0.4 





50 100 150 



0.8 



0.9 1 

Q./V 



a = 0.0 
a = 0.7 
a = 0.8 



1.1 



FIG. 4: The differential conductance as function of bias volt- 
age at hz — 10~* for different values of B2 and -B4. The 
easy-axis anisotropy of the molecule is D = 0.05 and the cou- 
pling to the leads J = 0.01. We have set Aq = 1. 



12D 




(+|5+|-) = (+|5+|->7^0 
(+|S-|-)^0 (+|5-|-)=0 

Jx Jy iJx Jy 5* 



FIG. 6: Derivative of the noise for the isotropic Kondo model 
without magnetic field a.t V = 50TR:,a=o. The asymmetry 
of the couplings is described by the asymmetry parameter a, 
where Jl,r = Jo(l±a). Inset: Noise for symmetric couplings 
and V = 50 Tk (solid line), V = 100 (dotted line), and 
V — 150 Ta' (dash-dotted line), showing a dip at S7 = itV. 

or B4, we get (+|S'+|-) ^ [see Fig. Mjo)], leading to 
JxJy > with a quantum phase transition at JxJy — 
to a Kondo effect The latter leads to an increased 
conductance at 1^ = ±h (see Fig. |4l). 



VI. NOISE FOR THE ISOTROPIC KONDO 
MODEL 



FIG. 5: Energy landscape and coupling of the states of a For h ^ T = and the isotropic case J' = J, we obtain 
molecular magnet with S = 7/2. The longitudinal anisotropy 

parameter D determines the height of the parabola —DS^. F = 7rJjjj|A=\/V, (81) 

The transverse anisotropy constants B2 and Bi couple every ^ 3r 

second or fourth state. For a pure B4, term, this leads to ^st ~ "T: 

(82) 

(+|S'+| — ) being zero or finite depending on the size of B4, 

since different states form the ground states I ±). — — Jj^^ f2|A=|y±f2| ± (83) 



with 51^, = 5ss', S^^J^ = Ss-s', and 

zs 

Vaa' ^ /-*" ^ f^a' i 

via''" = Ma - ^J■a' - s'h. (80) 

The decay rate and the stationary probability follow from 

r^Wn + w^iT> pf = Wn/r, and pf = Wi^/r^ Cm- 

rent and noise are obtained from Eq. (|79p as described 
at the end of Sec. IIVI 

Fig. |4] shows the differential conductance as function 
of the voltage at finite magnetic field for different values 
of the B2- and i34- anisotropy constants of a molecular 
magnet. For small B4 and i?2 = 0, we find no Kondo 
effect for specific spin values S' = 3/2 -I- 2m, m = 0, 1, . . . . 
In this case, the transverse exchange couplings J^^y = 
J(+|S'j^i ± "SmI— ) have the property J^Jy <Q according 
to ) = [see Fig. [5ja)]. By increasing either B2 



with 

^p-^ = ^inax(|y±n)|,r)'^d'/nd (84) 

in leading order, where Jnd = Jlr = Jrl and Jd = 
Jll = Jrr for a symmetric coupling to the leads. 
Whereas the decay rate and the current are cut off by 
the voltage, the noise is cut off by ± 17| which can 
be tuned to by setting — ±V . As a result, the 
noise is sensitive to the cutoff F of the couplings at these 
points (see Fig. [6|). There is a simple interpretation of 
the shape of the noise: It can be interpreted in terms of 
a golden rule expression with two superimposed currents 
^ "^nd f2lA=l^^±s^l l^i^^l with renormalized couplings (see 
Fig. [7]). For < y, the sum of both currents would be 
independent of 51 for bare couplings, but this balanc- 
ing does not hold for renormalized couplings and a dip 
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o<n<v o<v <n 

FIG. 7: Left: Interpretation of the finite finite frequency 
noise. It corresponds to the sum of two currents with dif- 
ferent voltages |y ± f2| and different renormalized couplings. 
If both couplings were bare, i.e., not renormalized and there- 
fore equal, the sum would be independent of for < ^1 < V . 
Right: For V < il, the sum increases with fl even for bare 
couplings because the absolute value of the currents is impor- 
tant. 

evolves (see inset of Fig. [5]), whereas for ft > V the noise 
rises again. The effect of T is prominent in the vicin- 
ity of = ±V. In Fig. H the derivative dS{n)/dn is 
shown. The shoulders around = V show the logarith- 
mic scaling of the coupling. By tuning the asymmetry, 
the relaxation rate F can be tuned which results in differ- 
ent cutoff heights of the shoulders, leaving the remainder 
of the noise untouched. 



VII. CONCLUSION AND OUTLOOK 

In this work, we have discussed a fundamental model 
of dissipative quantum mechanics: a local quantum sys- 
tem at fixed particle number coupled via spin or orbital 
exchange to several electronic reservoirs (the generaliza- 
tion to bosonic reservoirs is straightforward and goes 
along similiar lines). We have proposed that a micro- 
scopic derivation of cutoff scales from decay rates should 
be based on a formulation in terms of the reduced den- 
sity matrix of the local system, since the decay rates oc- 
cur naturally by their definition, namely, as the negative 
imaginary part of the eigenvalues of the kernel determin- 
ing the time evolution of the reduced density matrix. We 
have shown that a complete description of decay rates in 
RG formalism is only possible if one considers the full 
Keldysh structure, since relaxation and dephasing essen- 
tially arise from diagrams connecting the upper with the 
lower part of the Keldysh contour. Therefore, project- 
ing the RG equation from the very beginning on only 
one part of the Keldysh contour, one cannot obtain a 
microscopic description of decay rates. Although energy 
broadening terms might still lead to a cutoff of the pro- 
jected RG flow in this case, the physics of relaxation and 
dephasing is not included. Therefore, our approach pro- 
vides a consistent nonequilibrium RG formulation which 
can identify the generation of energy broadening, relax- 
ation, and dephasing at the same time, together with 



their influence on the RG flow of the vertices. 

Within our formalism, a particular problem arises due 
to the existence of an eigenvector with zero eigenvalue 
(the stationary state, which is always present and unique, 
at least in the absence of symmetry-breaking). It is a 
nontrivial technical issue to show that this eigenvector 
does not induce a flow to strong coupling. We have 
achieved this for a generic quantum dot in the Coulomb 
blockade regime by analyzing the one-loop RG equations 
in leading order, including all kinds of boundary vertices 
determining the dot distribution, the current, and the 
noise in frequency (Laplace) space. From a pure physi- 
cal point of view, one does not expect that the presence 
of a stationary state is correlated to the presence of a 
strong coupling fixed point, since the former is generic 
and the latter model-specific. Therefore, we believe that 
our leading-order analysis will hold in all orders but a 
general technical proof for this is still lacking. Neverthe- 
less, we have demonstrated within a one-loop formulation 
that decay rates cut off the RG flow generically. The 
precise prefactor of the various decay rates cannot be 
determined by our analytic formulation so far, since cer- 
tain irrelevant contributions of higher-order terms have 
been included in the one-loop equations, arising from the 
e^*^°*-factors in the definition of the interaction picture 
[see, e.g., Eq. ([4T|) ]. These factors are essential to pro- 
vide a cutoff scale for the generic case, but lead also to 
the unwanted effect that certain combinations of decay 
rates occur which prohibit an unambigious determina- 
tion of the correct prefactor. Also here, further devel- 
opments of nonequilibrium RG are needed to provide 
generic schemes for problems where the prefactor of the 
decay rates plays an important role, e.g., for problems 
with several decay rates differing by many orders of mag- 
nitude. 

The formulation of this work is useful for problems 
which stay in the weak-coupling regime with decay rates 
of the same order of magnitude. In this case, the dif- 
ferent prefactors will only lead to very weak logarithmic 
corrections. We have obtained the physically very natu- 
ral result that decay rates are only generated when the 
band width reaches the cutoff scales, defined by voltage, 
temperature, frequencies, or dot excitations. However, 
considering the Kondo model, when all these cutoff scales 
are zero, the only energy scale left is the Kondo temper- 
ature Tk, and we enter the strong-coupling regime for 
A < Tk- An interesting issue for future research is the 
investigation of the influence of decay rates in the strong- 
coupling regime. How or whether decay rates will cut off 
the RG flow also in this case is an open question and 
has so far not been analyzed. Due to the generic pres- 
ence of decay rates, we expect them to be important not 
only for weak-coupling problems. Such developments are 
highly desirable since no numerical method exists so far 
which is capable of dealing with the strong-coupling limit 
of dissipative quantum systems in a nonequilibrium sta- 
tionary state. Benchmarks for special problems are start- 
ing to be developed within the scattering Bethe ansatz 
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technique^^ but analytical and numerical methods for 
generic problems are still missing. 



The RG equation for Lu is obtained from Fig. ^h) 
and reads 
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APPENDIX A: DERIVATION OF THE 
RENORMALIZATION GROUP EQUATIONS 

In this appendix, we derive the RG equations (I56I) and 
([57)1 from Fig. [51 As explained in detail in Ref. \M, the 

A dependence of the vertices G^\f^ and the dot Liouvil- 
lian Ld have to be defined in such a way that the total 
sum of all diagrams stays invariant. This means that the 
derivatives of these quantities with respect to A have to 
cancel the corresponding RG diagrams of Fig. O Using 
Eq. (15^ , this gives for the vertex RG from Fig. ^a) 



f-lPlP-2 (^V2Pl _ . jPl JP2 tP2 jPl . 



r'P2Pl /-<PlP2 . jP2 jPl jPl j: 

^2'l',t/2 ^12,-t/2 • >^2' + ^l'-"^l+'^; 



Here, the contribution from the upper limit of time in- 
tegration can be shown to cancel exactly with a corre- 
sponding contribution arising from higher order (due to 
certain correction terms from time ordering; see Ref. 1 
for further details), i.e., there is never any divergence. 
Using 



PWl ■ 'Jl'if-'lrj 



jPi jPi . _ . 7P1 7P1 



(Al) 



together with Eq. (H^ and omitting the term : J^^J^^^_ 
on both sides, we obtain 



dA 



I "72-, 2 



ffPlP2 (~lP2Pl 

dA 12,t/2 ""2'l',-t/2 



-P2- 



" /2'-(-,2 



dA 



riP2Pi fiPi 
' ^2'\'tl2 ^Vl 



P2 

12 ^12,-4/2 



Inserting Eq. (j42|) for the contraction and taking matrix 
elements with respect to the eigenvectors of gives the 
RG equation ([TO)) for the vertex. 



dA 



i-iY / dt p[p'^ 

Jo 



f^PlPl f-lP2P2 . tPI jPl jP2 jP2 

^11' .t/2^22' .-t/2 ■ -'l+-'l'--'2+"'2'- 

Using Eqs. and (glD, we get 



i dt P2P2 



~dA 



_(^PlP2 -,PlP2 \^PlPl (^P2P2 

^^\'l'-,2+'i+,'^'-> ll',t/2 ^22',-t/2- 

Again, inserting Eq. (|42p for the contraction and taking 
matrix elements with respect to the eigenvectors of Ldj 
we obtain the RG equation ([57)) for the dot Liouvillian. 



APPENDIX B: CURRENT AND NOISE 

In order to calculate the probabilities, the current, and 
the noise from Eqs. ([^ . ([?T|) . and l[55|). one needs RG 

equations for the kernels Sn, S^, and Eq^ . The per- 
turbation series of So contains terms of the following 
structure: 



AqGG . . . GBq, 



(Bl) 



with the interaction picture of the two boundary vertices 
and defined by 



n.t 



iQt lEt 

e e 



APP 



-iLj:)t 



UPP _ -iilt iEt iL-Dt npp' 

-"ii'n.t — ^ e e -On/n, 



(B2) 
(B3) 



with E = ui — uj'i + iiai — IJ'a[ ■ Initially, Aq and Sji are 
independent of Q and given by the vertex G, defined in 
Eq. ([55)1 . For the single current kernels 'S'^ ^^'^ C^nY^ 
we have three different types of terms corresponding to 
whether the current vertex lies at the boundaries or in 
the middle of a diagram. 



y7 



-7t 



Al,^G...GBn 
AG...GGlG...GBa 



AG...GB, 



If 
n ' 



1-1 



AnG... GB.^ 
AnG...GG\^G...GB 



.GB. 



(B4) 



(B5) 



Here, X = Xn=o and all current vertices A^^ , A^^ , B^^ , 
B^ , and Gq are initially identical to the frequency inde- 
pendent current vertex G^, defined in Eq. (|49p . However, 
the interaction picture of all these vertices is defined dif- 
ferently and, therefore, they are no longer identical after 
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renormalization. With E = lui — uj[ + fiai — ^J■a[^ the 
various interaction pictures are defined by 

4/7, pp' _ iUt iEt All,Pp'~iLot 



(B6) 

e --e-^e-"^B^7;^^ (B7) 

^^^at^^Et^^L^tQ'y^,PP^^-^L^t^ (B8) 

A[y^^'=e^'A[y/e-^^-\ (B9) 



T>I1,PP' _ -iflt iEt „iLut tdIi,Pp' 

^ii'n,t 

^7,pp' 



^11' 



/7,pp' _ i_Et iLot f}ll,pp' 



(BIO) 



We note that the two current vertices A^^ and B^^ ac- 
quire only an imphcit frequency dependence via renor- 
mahzation but not an expUcit one from the interaction 
picture [see Eqs. (|B19P and (|B20|) below]. 

The RG equations analogous to Eqs. ([55)1 and ((Ff|) fol- 
low from Fig. [2] by replacing the vertex G by boundary 
vertices at the appropriate places, 



dA 



dE 



7t 



= An X 



/7 



(Bll) 

A X (B12) 



if7 X B, (B13) 



dAn 



dA 
dBn 
dA 



dA 



= An-G, 
= G ■ Bfi , 

= Al^-G + A-Gl, 



dB. 



77 



dA 
dA 



dA 
dA 



Gl-B 



n ■ -On, 



c c 

■ U + U ■ Uq, 



(B14) 
(B15) 
(B16) 
(B17) 
(B18) 
(B19) 
(B20) 
(B21) 



where we have used the abbreviations 



1 

[iX-Y)lf)^^ = i-dT / dt,5^,sign(c.2) 



{ 



(^^-^2)^4' ^^^^^ 



{X X y),fe 



d , 

« / dt—( 







dA' 



(^12,t/2)y (^21,-t/2)j'=- (^23) 

Note that we have already used the replacement —pf^^ 
\-fu,2 56'Tsign(tj2) inEq. (|B22| . As shown in Sec. [IVl 



this is justified in leading order. We note that for the 
calculation of the current and noise, given by Eqs. ([5T|) 
and the boundary vertices Aq, A^"^, and A^^ are 

only needed by summing over the Keldysh indices, 

X = Y, Xr-P' for X^AnA'^. A'^ . (B24) 



We find that the property is preserved under RG 
and holds also for the vertices A^ and B^ and for the 
kernel So, which by using Eq. ([55)1 reads 

^(0|Xff; = for X = Aa,Ba,G, (B25) 
pp' 

(0|A: = for A:==LD,En- (B26) 



From these properties, we get directly 
(OK' = (OIG^-', 



pp' 



\G}i^^ 



pp' 



l'-'f2=0' 



(B27) 
(B28) 
(B29) 



since these quantities fulfill the same RG equation and 
have the same initial condition. 

We now show in leading order that all RG equations 
are cut off by the decay rate F. The proof used in Sec. IIVI 
to show that the vertex G is cut off by F can be applied 
in a similiar way to the RG equation (jB2ip for the vertex 
G'li- For the boundary vertices, this proof does not work 
since the two vertices on the rhs of the RG equations 
(jB15P - (jB20P are not equal and the interaction picture of 
the boundary vertices differs from the one of G due to the 
absence of either e'-^^* (for An, AI^ , and A^^) to the left 



or e 



°* (for B^ , and B^') to the right. However, 
we can make use of the property (|B25[) to show that the 
eigenvector with eigenvalue zero cannot influence the RG 
equations in leading order. To show this, we consider as 
an example an RG equation of the form 



dXPP 
dA 



= {G-xyp 

^ (C-PP . j^p'p')(i) + {GPP ■ XPf)'^'^\ (B30) 



where the two terms on the rhs of (|B30|) correspond to 
the two terms on the rhs of Eq. (|B22[) , and the interaction 
picture of X is defindcd according to the boundary vertex 
B [for terms contributing to the boundary vertex A, the 
analysis is even simpler, since only the form summed over 
the Keldysh indices is needed; see Eq. I|B24[) ] . Defining 

X = ^XPP\ Xl = ^pXPp\ 

pp' pp' 

Xn = P'^""' ' = E PP'^'"' ' 

pp' pp' 
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= (G-X)(i) + 


(G-X)(2), 


(B31) 


= {Gl ■ X)('^ - 


K (G-Xl)(2), 


(B32) 


= (G • Xn)^'^ - 


H iGR-xY'\ 


(B33) 


= (Gfl-Xi)« 


+ (Gi-Xj^)(2). 


(B34) 



and the same for X ^ G, we get 

dX 
dA 

dA 
dXii 
dA 
dXLR 
dA 

We now consider the contribution when the two eigen- 
vectors i = j = in Eq. (jB22p have zero eigenvalue, 
i.e., Xi = Xj = (note that the eigenvalue cannot 
lead to any cutoff since the factor e~*^°* to the right is 
missing in the definition of the interaction picture of X). 
Due to (0|G = {0\X = 0, we see that this term does not 
lead to any contribution in Eqs. (|B3ip - (|B33p . Only for 
(IB34[) can the case i — j = contribute and no cutoff 
from a decay rate occurs. However, since the vertices Gl, 
Gr, Xl, and Xr are cut off by the decay rate, this term 
leads only to a logarithmic correction to Xlr which is in- 
tegrable and subleading (note that frequencies enter the 
argument of the logarithm, making this term finite when 
ingrated over the frequencies). Therefore, we see that we 
do not have to consider the eigenvalue zero in leading 
order (we expect that a similiar proof holds in all orders 
but this cannot be seen from one- loop RG equations). 

A similiar analysis can be performed for all terms of 
the RG equations (|BT5)) - (|B201) . With this resu lt, we can 
replace Ld ^ = I-^d, •] on the rhs of Eq. (|B22)) and 
introduce an overall cutoff factor 9r, 

(X.F)^i^ = i^9rdT I dtS^,sign{uj2) 



/ yPiP2 T^P2P'i I 
|_^12,t/2 ^ 21',-t/2 ^ 



VP2P1 YP1P2 
^2V,t/2 ^ 12,-t/2 



(B35) 



For Eq. (|B23|) . we use the same but neglect the decay 
rates on the rhs because this does not cause any diver- 
gence (the resulting 5- and principal value integrals are 
convergent) , 



XxY 



^ ^12,t/2 21,-t/2 



(B36) 



In the next step, we show that above all cutoff scales, 
the current vertices preserve their initial form (j49p to- 
gether with Eqs. ([55)1 and To prove this, we evalu- 
ate (jB35|) for A larger than all cutoff scales and get 

{X ■ r)^i^ = i (Xff ^ Y^'f'^ - X^'f'^ rfa^P^) . (B37) 

Inserting this form in the RG equations (|B15P - (jB21[) . we 
find that all vertices are independent of il and we get 



A = B = G and A^-* = A^^ = B^^ = B^^ = G^. 
Furthermore, inserting the initial form (|^^ together with 

Eq. dSSl) for the vertex G^i?^' = -^a,-y-6c,'^^)p'dpp.GPP 
on the rhs of the RG equation (|B2ip . we find 

j(^7,Plpi 1 1 ( 

-J^ = ^2'^\Pi(^ai-t " ^a2-t)[Gil^^ jG^lf^] 

— P'l(^a2 7 ~ ^a['y)[G2\f^ ,G{2^^]^ . 

Using the initial form ([M)) of the vertex G^^, we find 
[G++,G— ] = and we get 

11' — - -A,^r,A , r^pipi r^PiPi 



dA 



2 

d 

dA 



where we have used the RG equation ((6T|) above all cutoff 
scales in the last line. This shows that the initial form 
(j49p of the current vertex is preserved in leading order. 
Therefore, for all vertices X = G, A^ , , A^ , B^ , we use 
the form 



Yl,PP _ ^7 



fCI „' vPP 



(B38) 



(B39) 



also below the cutoff scales, together with 

vPP' _ X vPP 
^ll'O ~ "pp'^ll'Q 

for al l X = G,A,B,A^,B',A^,B^. Inserting the form 
(jB38[) for the current vertices into the RG equations 
(jB17[) - (|B21[) and neglecting all terms on the rhs which 
do not preserve this form (and become zero above all 
cutoff scales), we find the same RG equations for the 
7-independent vertices Aq, J5q, JLq, Bq, Gn but with an 
addtional factor i appearing on the rhs of the RG equa- 
tions. 



dA 


- ^{A'n 


•G H 


- A-Gn}, 


(B40) 


dBh 
dA 






- G_o • B} , 


(B41) 


dA'n 
dA 




•G - 


F An-G-n}, 


(B42) 


dBh 
dA 




Bh- 


\- Gfi ■ Boj , 


(B43) 


dGn 
dA 




•G ^ 


- G • Go} . 


(B44) 



Usi n g Eq. (jB35|l in the RG equations (jBlSjl . (|B16p . and 
(|B40)) - (|B44)) . one can easily prove the following useful 
relationships between the vertices [note that the form 
i^j' = [i/D, •] implies L'S' - (ij^')^]: 



Go^ 



Gr 



An^ = Bn, 
A'n = Bfi, 



(Ah) 



Bh- 



(B45) 
(B46) 
(B47) 
(B48) 
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Furthermore, we get the following properties for the ma- 
trix representations of the vertices: 



= 0, 



(B49) 



p S 



-iXrl:Z!^):,,,,s^s. (B50) 



for all vertices X = G, A, B, A' , , A' , B' . 

The RG equation (|B44[) for the vertex Go becomes 
especially simple. If one uses Eqs. (IB35|) and (jB39p and 
Lp' = [^^D, •], one finds that the initial form 



■ givn 



(B51) 



is preserved under renormalization. Therefore, one can 
project this RG equation exactly on the upper Keldysh 
contour and obtain 



dgn I r , . 



(B52) 



where in analogy to Eqs. (jB35P and (jBSp . we have defined 



{x ■ y)ii' = i-OrOr J dt 5^^sigii{uj2) 



{xi2,t/2V2V-t/2 + 2^21', t/2 2/12 -4/2} (B53) 



and 



AM iEt iHr>t „ ^—iHot 

givnj — e e e givne , 



(B54) 



with E ^ u!i — uj[ + jiai — A*a'j • Evaluating Eq. (jB52p 
gives Eq. ([75)1 of Sec. lIVI We note that this projection on 
the upper Keldysh contour is not exactly possible for the 
boundary vertices since the interaction picture is defined 
differently. Whereas for the vertex we can use 



r++ — n 



a similiar equation does not hold for the boundary ver- 
tices. Using Eqs. (jB45P and (jBSip . we get 



givn^ — 5i'i,-o- 



(B55) 



Finally, we note that for not more than two reservoirs, 
the generation of double-current vertices must not be 
considered, since they can be shown to be irrelevant, i.e., 
they are not generated above all cutoff scales, at least 
if the trace over the dot states and the sum over the 
Keldysh indices are taken [which is the quantity entering 
the noise formula ([55)) ]. To prove this, we use the leading 
order form ([B37[) above all cutoff scales and get for the 
RG of double current vertices G'^'^ 



Inserting the form ([B38)) . ([B39| . and ([B5T)l for the cur- 
rent vertex, we get after some straightforward manipula- 
tions 



dGll ^"'"^ _ 1 



dA 4A 

pip'i 

-{Sai-y' - Sa2l'){^a2'r " ^a[-()} Tro (512521' + 5'21'3l2) • 

For two reservoirs, the rhs of this equation can easily seen 
to be zero either for 7 = 7' or 7 7^ 7'. 

The complicated set of RG equations simplifies con- 
siderably if one considers a problem where the dot dis- 
tribution is diagonal {pn)ss' = Sss'{pn)s and where the 
dot eigenstates \s) do not renormalize (however, the dot 
energies Es can renormalize). This is, e.g., the case 
for the fully anisotropic Kondo model under consider- 
ation in this work, given by Eqs. ([21)) - ([26)1 . due to ro- 
tational invariance around the z axis. In this case, 
we need for the dot distribution ([^5)l . the current ()5ip . 
and the noise (|55p only the matrix elements X^s^s's' 
{X = Eq, Eq, Eq^, Eq^ ) for the kernels. As a conse- 
quence, we see from (|Blip - (|B14l) and (|B38p that only 
the components Xss.-- {X — Ao,^^,^^) and X..^ss 
{X — B^i, Bf^, Bf^) of the boundary vertices are needed. 
Using the matrix representation of Lq' = [^d,-] [see 
Eq. 



(-^D )sis[,S2s'2 — (£^si — Es'JSs^ 



S2 sL 



we find (Lg 



rcl\ 



{An)ss,- 



- and we get 
= (^r2)ss,-- = (Go)s 



(B56) 



(B57) 
(B58) 



since the interaction pictures and the RG equations are 
the same for the various quantities. For the boundary 
vertices A^ and Sji , we get this property only after sum- 
ming over the states and the Keldysh indices, since oth- 
erwise the second term on the rhs of Eqs. (|B42p and 
(|B43p contributes and leads to different renormalizations. 
Analogous to Eq. ()B29p . we get 

E(in'""')...- = E(G"o=lo)-,-- (B59) 

E(^/,'''^').... = E(GeLo)-...- (B60) 



Therefore, we need only the RG equation (|B52p for the 
vertex and we can easily evaluate Eqs. ()Blip - (|B14p 
by using Eqs. ([B36)) . ([B38)) . ([BST)) . and ([TO)) . 



PlP'l V\V2V\V'2 



|Lt^2 "-^21' ~ '-^21' "-^12 J • 



rf(Sa)ss,s's' ■ d in a \f- f+ 

dA --^rfA(^-^"')/-^' 

X < _ (Q^-nA, (B6I) 

n + LU-Lu' + y + ii] ^ J 
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for s ^ s' , and 



E 

s 

S 

E 

s 

^E 

s 

E 



dA 

[_ J7 + — + ?; + ir/ 

'^(^Sl''' )ss,s's' 



J J U) J ijj' 



(B62) 



* 2 "^la' "^la' (^i^ ^t^' ) /J"' 



+ — w' + y + i?] 



E' 



dA 



(B63) 
(B64) 



with y ^ lia - ^J.a' + Es ~Es^ For f7 = 0, Eq. (|B6T1) 
leads to Eq. ^ and Eq. (|Bg^ to Eq. ^ of Sec. [iVl 
giving the stationary dot distribution and the station- 
ary current [the time-dependence of the dot distribution 
and the current for an arbitrary initial state can also be 
calculated from Eqs. (|43)) and ([ST]) ]. The calculation for 
the noise simplifies considerably if the above matrix ele- 
ments do not depend on s (which is the case, e.g., for the 
isotropic Kondo model in the absence of a magnetic field 
due to spin symmetry). In this case, we can average over 
s and get from (|B64p and (|B26p the explicit formulas 



iEn^)ss,s's' = =0, (B65) 

S 

s 

= JE(^o=o)--'«'- (B66) 



Therefore, we get (nf2)ss,s's' = from Eq. ([43)) . and 
using Eqs. (|27p and ([55]) . we obtain for the diagonal noise 
[the nondiagonal one follows from Eq. ([29| ] for fl 



- ^E(^o" + siT,). 



(B67) 



1 i 

Z2n 



E(^o-s] 



E(^o=o) 



SS^S S ( J 



where we have used pf — 1 / Z with Z denoting the num- 
ber of dot states. Using Eq. (|B63[) . we get directly the 
RG equation ^ of Sec. [IV] for Y^ss'^^n + '>^T'n)ss,s' s' ■ 
Concerning the second term on the rhs of Eq. (jB67p . we 
use Eq. (|B62p and interchange <-> /^', w <-^- w' and 
s s' in the second term on the rhs of this equation. 
With the help of Eq. (|B55|) . this gives the result 



xP 



1 



r2 + W — + /ia — /io 



; j ^\g^,^,'n{uJ,u}')ss'\'^ 



Performing the integrals over uj and u' by neglecting the 
frequency-dependence of g^^iii{ijj,uj'), one finds that this 
term leads to a contribution of the order F x 0(4-1 4^)- 

^ Ao ' Ao ^ 

Therefore, it is irrelevant and is left out within our lead- 
ing order analysis. 
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